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Abstract 

We present experimental evidence of coexisting periodic attractors in 
a semiconductor laser subject to external optical injection. The coexisting 
attractors appear after the semiconductor laser has undergone a Hopf bi- 
furcation from the locked steady state. We consider the single mode rate 
equations and derive a third order differential equation for the phase of the 
laser field. We then analyze the bifurcation diagram of the time periodic 
states in terms of the frequency detuning and the injection rate and show the 
existence of multiple periodic attractors. 
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1 Introduction 



Semiconductor lasers have a wide range of applications because they are of rela- 
tively small size, they can be massively produced at low cost, and they are easy 
to operate. Applications of semiconductor lasers appear in many areas such as 
optical communication and high speed modulation and detection. Despite their 
successful technology, semiconductor lasers are quite sensitive to any external 
perturbation which may destabilize their normal output. A small amount of opti- 
cal feedback resulting from the reflection from an optical disk or from the end of 
an optical fiber is sufficient to generate pulsating instabilities. These oscillations 
are typically accompanied by higher intensity of frequency noise and affect the 
normal efficiency of the laser. 

Systematic experimental studies of semiconductor lasers, in particular long 
time series analysis of the intensity are not available because the time scale of 
the intensity pulsations is typically in the picosecond regime. Fourier spectra 
measurements show a gradual increase of oscillatory instabilities as parameters 
are changed but do not reveal what are the bifurcation mechanisms. Most of the 
progress on our understanding of these bifurcations comes from extensive numer- 
ical studies of simple models and their comparison to the experimentally obtained 
Fourier spectra. Specifically, these are the Lang and Kobayashi equations \^ for 
a single mode laser subject to optical feedback [@], the equations for two coupled 
lasers [Q], and the equations for a laser subject to external optical injection [Q]. 
The latter is presumably the simplest system which allows a study of instabilities 
induced by an external perturbation. 

The rate equations in normalized form for this system consist of two equations 
for the complex electrical field E and the excess carrier number N given by: 

E' = {l-ib)NE + r]Ei, (1) 

TN' = P - N - P{1 + 2N)\E\'^ (2) 

where prime means differentiation with respect to time t measured in units of the 
photon lifetime r^. The term r]Ei = t] ex.p{—i^lt) models the electrical field of 
the injected signal which is controlled by changing either its amplitude i] or its 
detuning ^l. The equations are the same for a laser subject to optical feedback or 
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for each laser in arrays of coupled lasers. They differ only by the last term in Eq. 
(|1]). The fixed parameters b, T and P are defined as follows: b is the linewidth 
enhancement factor (6 = 3 — 6) which measures the amount of amplitude-phase 
coupling; T = r^/rp is the ratio of the carrier and photon lifetimes (T ~ 10'^) and 
P is the dimensionless pumping current above threshold (|P| < 1). A series of 
recent experiments have explored a large variety of dynamical instabilities at high 
injection levels and have been simulated successfully by using Eqs. (|l]) and @. 
The comparison is particularly remarkable knowing the high simplicity of these 
equations. They correctly describe the observed period doubling cascade to chaos 
[|]], the period 2 bubbles as well as the details of the experimental map of all 
the instability regions in the detuning versus injection amplitude parameter space 

Of particular interest are the conditions for locking which corresponds to a 
steady phase of the laser field and occurs if the detuning is sufficiently low. This 
phenomenon is best analyzed by introducing E = _Rexp[2($ — fit)] into Eqs. ((ij) 
and @ and by studying the conditions for a steady phase $. If the injection level 
is very weak, it is reasonable to assume that the laser amplitude and the carrier 
number will approach their steady state values in the absence of injection (i.e., 
i? ^ 1 and ~ 0) and that the long time behavior of the laser is first described 
by the phase $. From the equations for R, $ and A^, it is then straightforward to 
find that $ satisfies Adler equation 

$' = fi-r/sin($). (3) 

Locking now implies the condition |fi| < t]. In addition to locking, Eq. (^ has 
been used to investigate the case of an unbounded phase (i.e., if |fi| >ri). Specif- 
ically, if 7] << fi, $ ~ fit which implies a small amplitude periodic modulation 
of the laser amplitude equation i.e., R' = RN + rj cos($) ~ RN + rj cos(fit). 

This case has been studied analytically [jSj] and is called four wave mixing 
because of the typical Fourier spectra. However, the practical interest of (^ is 
limited by the fact that T is a large parameter. The effect of T large can be first 
analyzed from Eqs. ([T]) and (^ with r/ = 0. We find that small perturbations from 
the steady state {R,N) = (1,0) slowly decays on an 0{T) time interval with 
rapid oscillations characterized by period 2n / lu = 0{\/T) where lu = ^J2PJT 
is the laser relaxation oscillations frequency. We take this property into account 
by introducing the new time s = cut and by reformulating Eqs. ([1]) and (@) in 
terms of normalized deviations from the steady state, namely, x = Nuj^^ and 
?/ = i? — 1 . By then investigating the limit — > (equivalently, T large) of the 
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resulting equations for x, y and $, we reobtain (Q) provided that = 0{r]) and 



r]<uj' = 0{T-^). (4) 

Thus, @ is only valid for relatively low injection. It explains the locking phe- 
nomenon and the entrainement of the laser in the four wave mixing regime but 
cannot lead to the multiple instabilities observed experimentally and numerically. 

In particular, we show in Figures la and lb experimentally obtained spectra of 
a semiconductor laser subject to injection operating at very low detunings between 
the frequencies of the slave and master laser. The spectra were taken at constant 
injection and detuning. They show the characteristic sidebands at the relaxation 
oscillations frequency indicating that a Hopf bifurcation has taken place from the 
locked steady state. The laser exhibited the two spectra at various times. The 
transition from one spectrum to another was abrupt strongly suggesting that the 
laser was visiting two different periodic states assisted by the intrinsic noise in the 
system. This implies that another stable periodic state may coexist with the Hopf 
bifurcation periodic state for the same values of the parameters. Figures 2a and 2b 
show spectra obtained from solving numerically the single mode rate equations ([l|) 
and (|^). The parameters used for these simulations were T = 1000, b = A, P = 
1.0, r] = 0.002 and il/cu = 0.1. These values of the parameters are of the same 
order as the values that have been used before by a number of authors to model the 
HLP-1400 semiconductor lasers (emitting at 830 nm). These lasers were also used 
in this set of experiments [Q]. Gain saturation, known to shift bifurcation points, 
did not alter our conclusions, producing only minor modifications of the spectra. 
Figure 3 shows the bifurcation diagram of the periodic states. In addition to the 
periodic states emerging from the Hopf bifurcation, there is an isolated branch 
of periodic states emerging from a limit point. The spectra shown in Figures 2a 
and 2b were computed at the points denoted by A and B in Figure 3, respectively. 
They compare well with the experimental spectra shown in Figures la and lb. 

Bifurcation diagrams exhibiting coexisting bifurcation and isolated branches 
of periodic attractors are not common for nonlinear bifurcation problems. Similar 
diagrams have been found numerically for the laser subject to optical feedback 



[|T7p suggesting that the emergence of coexisting branches of periodic attractors 
is a general feature of semiconductor laser instabilities. This motivates analytical 
studies of the laser rate equations which we describe in the next two sections. In 
Section 2, we formulate a third order phase equation which is a major simplifi- 
cation of the original laser equations. In Section 3, we analyze this equation and 
determine approximations for the steady and time-periodic solutions. 
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2 Formulation 



Intensive numerical studies of the first period doubling bifurcation in the simplest 
case of zero detuning [ [T0| ] showed that the linewidth enhancement factor plays a 
key role in the problem. Its relatively large value compared to the other constants 
in the problem explains why the intensity oscillations are nearly harmonic in time 
and how a strong coupling between phase and laser amplitude is possible. This 
can be demonstrated by a new asymptotic analysis of the original equations us- 
ing b as the order parameter. To this end, we first introduce the decomposition 
E = i?exp[i($ — ^It)] into Eqs. (|T]) and (0) and reformulate the resulting three 
equations in terms of new variables defined by 

s = ujt,x = Nbcj-^andy = b(R - 1). (5) 

In (IJ), time is scaled by the relaxation oscillations frequency u = a/2P/T and 
X, y represent deviations from the zero injection steady state solution (R, N) = 
(1, 0). The resulting equations for x, y and $ are then analyzed for large b assum- 
ing uj = 0{b'^) (see Appendix). The reduced problem - accurate to 0{b^^) - is 
described in terms of the phase $ and is given by (p5|), or equivalently, 



+ - 2b-^A)<i>" + $' - A - A cos($) 
= [A sin($)$'2 _ ^ gjj^($) _ 2$'$"] W 

where prime now means differentiation with respect to time s. The parameters A, 
A and ^ are defined by 

rib ^ VL , 1 + 2P 

A=-,A = - and ^ = oo—— (7) 
UJ oj 2P 

and correspond to the normalized injection, normalized detuning and laser damp- 
ing parameters, respectively. The leading order problem as 6^^ ^ (tu small but 
fixed) is 



+ ^<l>" + $ - A - A cos($) = (8) 

and was previously examined for the case A = (zero detuning) [pDj]. In par- 
ticular, its validity has been carefully investigated by comparing the numerical 
bifurcation diagrams obtained from Eqs. ((T]) and @) and from Eq. (§). The two 
bifurcation diagrams have been compared for various ranges of values of the pa- 
rameters T and b including the experimental parameters estimated in [|]], [§] and 
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[{7p. The agreement between exact and approximate diagrams is excellent pro- 
vided that the injection level is not too large (i.e., A must be 0(1) or smaller). For 
the case A 7^ 0, a direct comparison between experimental data and solutions of 
Eq. (|]) has been undertaken if A is close to 2 and in the four wave mixing re- 
gion. If A ~ 2, the laser exhibits two period doubling bifurcations which are not 
observed for other values of A and motivated the interest for this particular case. 
The optical spectra have been obtained for A below and above 2 which allow us 
to study both the sudden transition to subharmonic resonance and the progressive 
shift of the frequencies [[TT]]. 

In this paper, we investigate Eq. (^. The steady states and their linear stability 
properties are easily determined for small b^^ and cu = 0{b^^). We find that the 
steady states emerge from limit points located at 

A = A (9) 

and that they may change stability at Hopf bifurcation points located at 

A = yA2 + (^-26-iA)2 + 0(6-2). (10) 

From @ and {^^, we note that a limit point and a Hopf bifurcation point may 
coalesce if A = A^ and A = A^ where 

Ac ~ ^&/2andAc = Ac. (11) 

This degenerate Hopf bifurcation point was already noticed from the exact linear 
stability analysis [|12p. It corresponds to a zero eigenvalue and a pair of imaginary 
eigenvalues of the linearized theory and may lead to a secondary bifurcation to 
quasiperiodic oscillations [|13]]. 

3 Bifurcation equations 

The phase equation (^ is a major simplification of the original laser equations. It 
models the laser problem as a weakly damped harmonic oscillator strongly driven 
by the phase of the laser field but it is still too complicated for exact solutions. In 
this section, we determine an asymptotic solution of Eq. (|^) valid for small b~^, 
A, A and ^ by the method of multiple scales [[T5|]. To this end, we expand the 
parameters as 
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A = 6-i(A + b-'Xi + ...), A = h-'{6 + b-i<5i + ...), 



and seek a solution of the form 



<l>(s, r, 6-1) = $0(5, r) + 6-i<l>i(s, r) + ... (13) 

where r = fe^^s is defined as a slow time variable. The fact that the solution 
depends on two independent time variables implies the chain rule $' = ^g+b^^^r 
where subscripts mean partial derivatives. 

Inserting ( [T^ and ( [131 ) into Eq. (^ and equating to zero the coefficients of 
each power of leads to a sequence of linear problems for the coefficients 
$0; ••• Solving the equation for $0 gives 

= Asm{s + v) + B + 0{b-^) (14) 

where the amplitudes A and B and the phase v are functions of the slow time r. 
We obtain equations for A, B and v by considering the problem for $1 and by 
applying solvability conditions [jT^] . The equation for v is simply v' = which 
means that v is equal to its initial value. The equations for A and B are given by 

A' = --aA + Xsm{B)Ji{A) (15) 
2 



B' = 5 + Xcos{B)Jo{A) (16) 

where Jo{A) and Ji{A) are Bessel functions which come from the expansion of 
cos($) in Fourier series. Note that the modified phase equation (^ and the leading 
order phase equation (||) lead to the same amplitude equations at this order of 
the analysis. Eqs.([T3|) and ( [T^ ) admit steady state solutions which we analyze in 
detail. The simplest solution is given by 

(1)A = and A = -A cos -^(B) > (17) 

and corresponds to the steady state of the original equation (||). The second steady 
state solution of Eqs.(P3|) and ( [T^ is characterized by the fact that A ^ and 
corresponds to time-periodic solutions of Eq. (|8[). It is instructive to examine this 
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solution for A = first. From Eq. (|T^ with B' = S = 0,wo find two possibilities. 
The first family of solutions verifies the condition cos(i?) = and is described by 



TT 



(2a)B = ±- and A 



"2Ji(A) 



> 0. 



(18) 



From ([T8|), we determine a Hopf bifurcation branch emerging from Solution (1) 
at (A, A, B) = (A^, 0, f ). A^ is found from dH) by taking the limit A ^ 0. We 
obtain 



A 



H 



(19) 



which is correctly matching ( [T0| ) if A = 0. But Eq. ( [181 ) describes other branches 
of solutions. These branches of solutions are isolated and emerge from limit points 
that satisfy the condition A' (A) = 0. Using the second expression in ([T8|), we find 
that these limit points are located at (A, A, B) = (A^, A^, ±|) where A^ and A^ 
satisfies the condition 



MA' 



and A^ 



^A' 



"2Ji(A^) 



> 0. 



(20) 



The second family of solutions characterized by an A 7^ verifies the condi- 
tion Jo{A) = and is described by 



(26) Jo (A) = and A 



2sin(B)Ji(A) 



> 0. 



(21) 



This solution emerges from Solution (2a) at the pitchfork bifurcation point located 



at (A, A, B) = (A^^, A^^, ±f ) where 



Jo(A^^) = and A^^ = ± 



2Ji{APF) 



> 0. 



(22) 



In the case of nonzero detuning (A 7^ 0), we find that all A 7^ steady state 
solutions of Eqs. ( [151 ) and ( [T6| ) are described by (in implicit form) 



A 



A 



Jo (A) 



^A 



2Ji(A) 



1/2 



and B = arccos 



A 



A Jo (A) 



(23) 



Eq. ( P3| ) reveals several branches of solutions which converge to Solution (2a) 
and Solution (26) as A ^ . In Figure 4, we show the bifurcation diagram 
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of the solutions of Eqs. (|T5]) and (|T^ in terms of A vs A. The solid and dotted 
lines correspond to A = and A = 0.01, respectively. All stable and unstable 
branches of solutions have been drawn. Note the coexistence of multiple branches 
of solutions emerging from limit points. If A = 0, these limits points correspond 
to simple roots of Bessel functions (see Eq.(pO|) and Eq.(|2^)). 

The numerical bifurcation diagram of the stable solutions of Eq.(H) with A = 

is shown in Figure 5. Figures 5a and 5b represent the extrema of the phase $ 
and the intensity variable y = b{R — 1), respectively. The pitchfork bifurcation at 
( ^2[ ) leading to two distinct branches of B is clearly identified in Figure 5a. These 
two branches are not seen in Figure 5b because ?/ = $" ~ —A sin(s + v) from 
( p3[ ) and is independent of _B, in first approximation (equivalently, the extrema of 
y are equal to ± A, in first approximation). The figure shows one isolated branch 
of stable periodic states that coexist with the bifurcation branch. The bifurcation 
at A ~ 0.65 is a period doubling bifurcation which is not predicted by Eqs. ( [T3| ) 
and ( [T^ and which requires a different analysis [|TOp. 

The unfolding of the A = bifurcation diagram as A is progressively in- 
creased is illustrated by the dotted lines in Figure 4. In Figure 6, we represent 
the extrema of y for all stable solutions determined numerically from (||) with 
A = 0.01. The agreement between numerical and analytical solutions is excellent 
(compare Figure 4 and Figure 6). 

The amplitude equations (|T5]) and (|T^ admit a third family of solution if 

1 A/ A I > 1. These solutions correspond to A = and an unbounded B satis- 
fying Adler's equation, B' = A + Acos(-B). We have verified that all solutions 
of this equations are linearly stable with respect to small perturbations in A. $ 
is then unbounded but $' is bounded (equivalently, a; ~ A — $' and y ~ $" are 
bounded). This solution is the four wave solution for small A and A = 0(A). 

In summary, our analysis of the phase equation (||) has revealed several branches 
of periodic solutions which allows the coexistence of different periodic states os- 
cillating harmonically in time with similar periods but exhibiting differing ampli- 
tudes. Coexisting periodic states oscillating harmonically in time were also ob- 
served for the laser subject to optical feedback [[T7|]. Both for the laser subject to 
injection and the laser subject to optical feedback, it is the strong phase/amplitude 
coupling parameterized by the linewidth enhancement factor b that is the main 
mechanism leading to this multiplicity of attractors. 

Note that these periodic states appear as isolated branches of solutions. The 
traditional way to determine the bifurcation diagram of the long time regimes is 
to follow each branch of solutions as a parameter is progressively changed and 
carefully determine changes of stabilities between each states. This method is 
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used both experimentally and numerically but will fail to find isolated branches of 
solutions. 

Isolated branches of solutions have been found numerically and experimen- 
tally for periodically modulated lasers [ |T^ ] . These branches correspond to differ- 
ent resonances and coexisting periodic states admit different periods. In the case 
of a semiconductor laser subject to feedback, two periodic coexisting attractors 
were observed experimentally and computed numerically by Mork etal^. More 
recently, De Jagher et al [ p^ reported two coexististing periodic solutions which 
may result from the emergence of isolated branches of solutions as described in 
this paper. From a practical point of view, it could be interesting to force the semi- 
conductor laser to operate at these higher intensity states and investigate possible 
application in communications or in logic gates. 

4 Discussion 

We have shown analytically that the single mode equations admit multiple peri- 
odic states emerging from limit points in addition to periodic solutions that appear 
as the result of a Hopf bifurcation from the steady state. In general, the spectra of 
two of the limit cycles might appear very similar if they are taken at an injection 
level at which the amplitude saturates and approaches constant values. However 
for low injection levels and away from saturation the spectra can exhibit dramatic 
differences. This is clearly demonstrated by the computed spectra shown in Figure 
2. In particular small perturbations in the phase of the field can induce dramatic 
changes in the optical spectra. This is precisely the reason that the multiple insta- 
bilities exhibited by the phase of the field can be clearly observed in the spectra. 
Such effects include center line suppression, in which the center line of the laser 
is almost devoid of energy most of which is carried by the sidebands at the re- 
laxation frequency. This is a typical feature of a laser with a large value of the 
linewidth enhancement factor. 

The effect of gain saturation can also be incorporated in the context of the 
phase equation formalism. Its main contribution is to shift Hopf bifurcation and 
limit points but it would not modify the qualitative properties of the bifurcation 
diagram. 

We note that the period doubling bifurcation point can also be captured analyt- 
ically from Eq. (|^) by seeking a solution of the form $ ~ A sin(s) + D cos(s/2). 
If A = 0, we obtain 
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KpD ^ 0.62 (24) 

which compares well with the prediction of Figure 5a. If A 7^ 0, we find that the 
period doubling bifurcation point Ap/)(A) is - in parametric form - given by 

where O(^) terms have been neglected. 

The existence of sustained quasiperiodic oscillations has been suggested nu- 
merically but require a more detailed bifurcation analysis. Furthermore, experi- 
mental evidence of quasiperiodic oscillations is still lacking. This is partly be- 
cause the laser is not stable at this injection region and frequently mode hops to 
a different cavity line. In addition the transition to the quasiperiodic behavior is a 
weak phenomenon and it is difficult to extract the signal from background noise. 

In conclusion by taking advantage of the two large parameters T and h that 
are naturally present in semiconductor lasers, we are able to reduce the original 
single mode rate equations into a third order pendulum type phase equation. This 
equation then allowed us to predict isolated branching of solutions which can not 
be anticipated using traditional continuation methods. 



5 Appendix: derivation of the phase equation 

We first introduce the decomposition E = i?exp[z($ — Qt)] into Eqs. ([T]) and (|^) 
and obtain the following three equations for A, $ and 

R' = NR + r]cos{^), (26) 

= fi-6iV- ^sin($), (27) 
R 

TN' = P-N -P{1 + 2N)R^. (28) 

In Eqs. (|2^)-(|2^), T is a large parameter that multiplies A^'. We remove this 
source of singularity by introducing the new variables s, x and y defined by (||). 
Eqs. (|2|)-(|2|) then become 

y' = x{l + b-^y) + Acos{<l>), (29) 
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^' = A-x-Ab' 



_i sin($) 
1 + 



(30) 



x' = —C,x — y — 2ub ^xy — -b ^y"^ — ub '^xy^ . (31) 

where ^, A and A are dimensionless parameters defined by Note that ^ is 
proportional to uj and that uj is an 0(T^^/^) small quantity. Our problem depends 
on two small parameters, namely b^^ and uj. The derivation of our simplified 
problem will be based on the limit b^^ assuming uj = 0{b^^). 
^From Eqs. (P^), we find x as a function $ 



X = A - $' - Ar^ sin($) + 0(r2) 
and from Eq.(^T]), using (P2|), we find y as 

y = -^x-x' + O(r^) = $" + 0{b-\uj) 



(32) 



(33) 



(Note: the 0(6^^) correction term in ( P^ is needed but the 0(6^^) and 0{uj) 
correction terms in ( p3| ) are not). In order to derive an equation for $ only, we 
first differentiate Eq. (|3T|) once and eliminate y' using We find 



-^x' - [x(l + r V) + A cos(<l>)] - b-^y {x + A cos($)) + 0(6"^ 6"^ 



(34) 



We next use ( |32| ) and ( |33| ) and eliminate x and ?/ in (^. This leads to the follow- 
ing third order equation for $ given by 



$' - A - Acos(<l>) 



+6-1 [$" (6^ - 2(A - $')) + A sin($) - A sin($)$'2] + 0{b-^, b'^uj) 
where 6^ = 0{buj) is 0(1). 



(35) 
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Figure captions 



Figure 1 Experimentally obtained optical spectra of the laser subject to injection 
in the neighborhood of the Hopf bifurcation. /,. is defined as the free running 
relaxation oscillation frequency, fr = -^ \/^^PJT. 

Figure 2 Numerically computed optical spectra from the single mode rate equa- 
tions (|1|) and (0) for T = 1000, 6 = 4, P = 1.0, rj = 0.002 and n/uj = 0.1. The 
zero of the frequency scale denotes the frequency of the master laser. 

Figure 3 Bifurcation diagram of the extrema of the deviation of the field from 
steady state as a function of injection. The solutions are computed from the rate 
equations (|T]) and (Q) using the same values of the parameters as in Figure 2. The 
points denoted by A and B indicate where the spectra in Figures 2a and 2b were 
numerically computed. 

Figure 4 Bifurcation diagram of the amplitude A of the periodic solution versus 
A constructed form the analytic solutions of the phase equation for T = 1000, 
b=10,P = 0.5. Solid lines are for A = and dotted lines for A = 0.01. 

Figure 5 Bifurcation diagram obtained numerically from the phase equation (||) 
for T = 1000, 6 = 10, P = 0.5, and zero detuning, (a) shows the extrema of the 
phase $ versus injection A; (b). shows the extrema of the field versus injection. 

Figure 6 Bifurcation diagram obtained numerically from the phase equation @ 
for T = 1000, 6 = 10, P = 0.5, and A = 0.01. 
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